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ABSTRACT 

Understanding the progenitors of core collapse supernovae and their popu- 
lation statistics is a key ingredient for many current studies in astronomy but 
as yet this remains elusive. Using the MESA stellar evolution code we study the 
dependence of the lower mass limit for making core collapse supernovae (SNe) 
as function of initial stellar metallicity. We find that this mass limit is smallest 
at [Z] « -2 with a value of ~ 8.3 M . At [Z] = the limit is ~ 9.5 M and 
continues to rise with higher metallicity. As a consequence, for a fixed initial 
mass function the supernova rate may be 20 % to 25 % higher at \Z\ = —2 than 
at [Z] = 0. This affects the association of observed SN rates as a probe for 
the cosmological star formation rate, rate predictions for supernova surveys, and 
population synthesis studies. 

Subject headings: supernovae: general — stars: evolution — stars: formation 
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Introduction 



Most massive stars end their lives as core collapse supernovae (CCSN; IColgate fc White 



19661 ). There is a lower mass limit below which a single star of a given metallicity, Z, wi ll 
not undergo such a supernova (SN; lEldridge fc Toutl 12004 ; iHeger et al.ll2003l ; ISmarttl 120091 ) . 
Below this limit, a star ejects its envelope as an asymptotic giant bran ch (AGB) star, forms 
a planetary nebula, and becomes a white dwarf ( lEldridge fc Toutll2004l ); the remaining mass 
of the star will be less than the Chandrasekhar mass, Mch ~ 1.46 Mq. Above the limit, 
however, the star buil ds up an iron core by nuclear fusion until the mass of the core exceeds 
a critical mass, M crit (lEldridge fc Toutl 12004 ; IHeger et al.ll2003l ). In such a star core collapse 
occurs: the iron core will fall inwards until repulsive nuclear forces stop its collapse. 



The actual value of this lower mass limit depends greatly on a variety of factors (INomoto fc Hashimoto 



19881 ). including the numerical simulations code used and how the relevant physics is modeled, 
e.g., the mass loss rates, mixing processes, convective boundary layers and semiconvection, 
which change t he resulting interior co mposition structure of the star. Rotation may also 
affect the limit (lEldridge fc Toutl 120041 ). but be we do not explore this in the p resent study. 
Meta llicity- independent values of 8 M — 10 M Q are often used in the literature (IHeger et al. 



20031). though we expe ct that it should be affected by metallicity (jCassisi fc Castellani 



1993 



Eldridge fc Toutl 120041 ) . Current models (lEldridge fe Toutl 120041 ) of the lower mass limit as 
a function of Z have been made with only a limited resolution in Z and cover a rather 
limited Z regi me. Observat ional determinations are too uncertain to explore the lower 



limit precisely. ISmarttl (120091 ) determines from a combination of supernovae detections and 
the Salp eter initial mass fun c tion ( IMF) that the lower limit lies between 7M and 9M Q , 
whereas ICassisi fc Castellanil (119931 ) focused on the dependence over a wide range of Z val- 
ues. These investigations, however, have not provided detailed information on how the lower 
limit depends on Z. 

A star with initial mass lower than that requi red for a class ical CCSN may also un- 
dergo core collapse triggered by electron captu res dNomotol Il984 ); this is expected to oc 



cur for a narrow range of initial stellar masses (jMiyaji et al. 



1980 



Mivaii fc Nomotolll987 



Poelarends et al.ll2008l ). iPoelarends et all (120081 ) estimates this range to be between 0.25 M 
and 0.65 M Q . We, however, do not explore this case; we instead focus on the minimum mass 
required for classical core collapse events. We also do not consider binary stars in which 
the mass limits will be altered by interaction including mass transfer and accretion. The 
minimum mass required for a star to still be able to ignite carbon burning and possibly 
become an electron capture SN is denoted M up , while that required for a star to undergo a 
classical core collapse event not triggered by electron captures with an ONe degenerate core 
is denoted M up '. Thus, we focus on M up '(Z). 
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We present results from a grid of non-rotating stellar models that map the minimum 
supernova mass as a function of metallicity. In £J2] we describe the models used and the 
physics relevant to supernova progenitors. Then in §3] we give the results from the grid of 
stellar models and provide a metallicity-dependent fitting function for the minimum mass for 
supernovae, M up ' (Z). We discuss the implications of the results on cosmology and galactic 
chemical evolution due to changes in the number of SNe with M up '(Z) in §H The discussion 
includes a brief summary of uncertainties and physical effects that cause the observed trends. 
Our conclusions are given in £0 



2. Setup and Simulations 



We used MESA (Modules for Experiments in Stellar Astrophysics, http : //mesa, sourceforge .net/) 
code revision 3290. MESA is a modern, open source package for computational astrophysics 
(IPaxton et al.l 120111 1 and includes a ID stellar evolution code, MESA star whose parts have 
been tested internally and verified with well-known evolution results. 

We computed a grid of non-rotating stellar models, with compositio n scaled linearly ac- 
cordin g to Y = 0.24 + 2Z. For the rel ative mass fractions of metals we use lGrevesse fc Sauval 



(119981 ) and I Anders fc Grevessd fjl989l ) for Li, Be, and B. The computation grid covers 24 val- 
ues of Z between and 0.04, and the initial masses ranged from M = 8.2 M Q to M = 10 M Q . 
For each model we follow the evolution from the pre-main sequence to the point where we 
can determine whether a core collapse supernova should result. 



Whereas recent solar abundance determinations (lAsplund et al.l 120091 ; lLodders et al. 



20091 ) have different abundance ratios and, foremost, a different absolute metallicity, this 
should have little effect on our overall conclusions. Additionally, the use of scaled solar 
metallicity is a common simplification, and whereas it does not account for details of the 
galacto-chemical evolution of the different species, nor the spread of abundance ratios for 
a given metallicity for different environments or times, it should still approximate the key 
properties of how varying metallicity changes the supernova mass limit. 



2.1. Mass loss 



Mass loss can strongly affect stel lar evolution and the final fate of the star, especially for 
hi gh mass and high me tallicity (e.g.. lHeger et al.ll2003l ). Here we applied the mass loss rates 



of de Jaeer et al. ( 1988h to the ent ire range of stars. Despite their age they are still considered 
adequate (jEldridge fc Toutl 120041 ) . These rates are empirical but they agree well with the 
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theoretical rates of IVink et all ( I2000l . l200ll ; Eldridge fc Toutl 120041 ) . The mass loss rates 
were scaled by a constant efficiency factor rj = 0.8 to adjust the empiri cal rates such that 



evolu tionary tracks for non-rotating models match observational data ( IMaeder fc Meynet 



200l|). 



We scale the mass loss rates with metal licity using M (Z n ) = M ( Zp) \J Z /Z Q ( Kudritzki fc Puis 



20001 ) where we use the s ame = 0.019 ([Anders fc Grevessdll989l ) as was used in the un- 
derlying mass loss rate of de Jager et al. ( 19881 ). and Z n is the initial metallicity of the star. 
We added this scaling to MESA star for our study. 



2.2. Overshooting 



Convective overshooting increases the size of the helium core (jSchroder et al.l Il997f ) 
and thereby affects the mass limit for supe r novae . We use the exponential diffusive model 
based on the prescription of iFreytag et al.l (119961 ). The diffusive coefficient, Dov, is given 
by -Dov = Do exp{— 2z/ (fHp)}, where z, Hp, and Do are determined by the particular 
overshooting region: z is the location from its boundary, Hp is the pressure scale height at 
the convective boundary, and Lb = t^ ase t. Here f base is a typical velocity a t the boundary 
of the convective zone, t = Hp/v^e a typical timescale (IFreytag et al.lll996l ). and / is a free 
parameter that governs the efficiency of th e overshoot mixing. For all models we adopt a 
value of / = 0.016, following iHerwigl ( 12000l ). 



2.3. Determining the fate of the star 



The evolution past central helium burning of a star with initial mass between 8 M and 
10 M generally proceeds along the AGB. This b egins by dredging up the helium sh ell above 
the core before usual helium shell flashes set in (jHerwigil2005l ; llben fc Renzinilll983l ). A star 
that produces a CCSN, however, will be halted in this evolution usually before the shell 
flashes by igniting neon burning. Once neon has ignited, the core typically burns all the way 
to iron and a core collapse ensues. 



The maximum mass of a stable white dwarf is Mch = 5.83 /i c 2 M (j Chandr asekhar 

1935h . If the core consists entirely of 12 C and 16 0, then /x e = 2 and M Ch = 1.46 M . But 
if the core is composed of 56 Fe, then /x e = 2.15 and M h = 1.26 M Q . Hence it is necessary 
to consider the composition of the core to obtain a useful value for Mch- For our models we 
have found that a value M h = 1.38 M was a good indicator for ig niting advanced burnin g 
phases, and this value is consistent with those used by other studies (jEldridge &: Toutll2004l ). 
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It corresponds to a /i c ~ 2.055. 



Eldridge &: Toutl (120041 ) find that the final fate of the star depends greatly on the second 



dredge-up, since the second dredge-up directly impacts the mass and chemical composition 
of the core by removing material of a certain composition from the core. The interplay of 
dredge-up and burning processes determines whether the cores can reach Mch and hence the 
final fate of the stars. For example, for stars of about 7M — 9M Q initial mass, dredge-up 
reduces the core size below Mch so that a SN will not result; for stars of initial mass > 11 M 
dredge-up will not occu r, and nuclear burning i ncreases the core mass until it cannot support 
itself and a SN results feldridge fc Toutl l2004h . 



Therefore simulations were run at least until the second dredge-up was finished by the 
time we determined the final fate. We assumed the outcome will be a white dwarf (WD) 
when: i) e u > e nuc in the core, i.e., energy loss due to neutrino emissions is greater than energy 
generation by nuclear burning, so that the core continues to cool; and %%) M corc < 1.38 M Q . 

These two criteria imply that further burning in the core will not ignite. Here we do not 
consi der the case of electron capture supernovae (ECSN) in super- AGB stars (IPoelarends et al 
20081 ). In that case, an O Ne core collapses due to electron captures o n 20 Ne and 24 Mg be- 



fore onset of Ne burning (jMiyaji et al. 



19801 : Mivaii fc Nomotolll987h . Those would allow 



supernovae from stars with initial masses below the lower limit for core collapse supernovae 
as studied here. 

Conversely, we use the following criteria as indicators that a star becomes a core collapse 
SN: i) The onset of Si burning, the final stage before core collapse; and ii) M core > 1.38 M Q . 

It may not be necessary to wait for Si burning, since after the onset Ne burning the core 
should proceed to Si burning. To assure this, our model s were run to Si burning , nevertheless. 
Since the Si burning timescale ranges from about 18 d flWoosley fc Jankall2005l ) to about 1 yr 
and since this is much shorter than the dredge-up timescales, the core mass will not shrink 
further. The second criterion ensures that core is large enough to undergo core collapse. We 
require that the stellar core must be well established and therefore the models must be run 
at least until they finish the second dredge-up. 



We stopped the evolution simulation when both criteria for either one of the stellar fates 
(WD or SN Type II) were met and allowed us to determine the fate of the model. 



-6 - 



3. Results 

In Table [H we list for each metallicity, Z, the minimum mass M up ' of a model that 
produced an SN Type II. A metallicity-dependent transition mass, M u p (Z), can be fitted 
by 

M u P ' (z) f Elo^ :[2]>-8.3 

M 

with coefficients aj given by 



±0.15 



9.19 : [Z] < -8.3 



a = +9.40858 
ai = +1.14548 
a 2 = +3.96411X10- 1 
a 3 = +2.96185xl0~ 2 



a 4 = -8.79237xl0~ 3 

a 5 = -1.96134xl0~ 3 (2) 
a 6 = -1.11999X10 -4 . 



For [Z] > —8.3, where [Z] = logZ/Z Q , this formula has a 1 a deviation of only 0.056 M Q 
(0.112 for the entire range). 

We find that for [Z] < —8.3 the lower mass limit becomes metallicity-independent and is 
best fit by a constant value of 9.19 M Q . Figure [TJ shows the fate for each model we computed 
as a function of its initial M and Z and our fitting function M np ' (Z). 



4. Discussion 



Our results are con s istent with both ob servations ( Smarttl2009 ) and the current paradigm 
( jEldridge &: Toutl l2004t iHeger et al.l 120031 ) that the minimum supernova mass lies between 
8 and 10 M Q . Specifically, previous findings (jEldridge &: Toutll2004l ; iNomoto fc Hashimoto 
19881 ) that the fate of the star depends on the evolution during the second dredge-up are con- 
firmed. For all models, an ONe core grew during dredge-u ps. We could determina te whether 
the star makes a supernova after the second dredge-up. lEldridge fc Toutl ( 120041 ) note that 
this behavior is because of early C burning, which is confirmed by our models, including 
off-center ignition of C burning. In a few cases the more advanced burning stages (usually 
Si burning) did not occur. In those cases we assumed that the star ejects its envelope and 
peacefully becomes a WD. Once Si burning ignites, however, the star will evolve to core 
collapse. 
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Fig. 1. — Fate of stars as a function of mass and metallicity. Filled green circles indicate 
models that make core collapse supernovae, red crosses show models that do not, and the 
black curve indicates our fit. 
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4.1. Uncertainties 



We do not find a sharp transition mass in our simulations. Instead, for a given Z, we 
find a switching back and forth between making SNe and WDs several times as the initial 
mass is increased. Most of this might be attributed to numerical noise. We note, however, 
that the non-linear nature of stellar structure, e.g., the interaction of different shell burning 
phases during carbon burning, may cause similar phenomena. But we cannot quantify the 
relative magnitude of this effect, if present. It is included in our error estimate. In "real" 
stars, additionally, convection will be chaotic-there is "weather" inside stars-but this is an 
effect beyond the realm of our simulations, next to rotation, binary stars, etc. Based on our 
detailed studies, e.g., for [Z] = —3.7, we estimate an uncertainty of the transition mass of 
about ±0.15 M for all values of Z, excluding systematic errors due to the model. 



4.2. Impact of Metallicity on Evolution 

At the highest metallicities we see the effect of increasing mass loss (§ 12. ip that reduces 
the mass of the star and hence shifts the mass limit upward. Additionally, an increase of 
opacity due to higher metal contact can also lead to more efficient dredge-ups that reduce the 
size of the core. The increase of the transition mass at low metallicities ([Z] < —2), where 
mass loss becomes unimportant, is dominated by a decrease in opacity. Therefore a smaller 
fraction of the core will be convective, leading to larger transition mass as [Z] decreases below 
—2. Below [Z] < —8 the initial abundance of CNO isotopes is so low that it is insufficient to 
efficiently drive hydrogen burning: the star contracts until carbon is made by fusing helium 
in a primary way, then re-expands, hence the evolution becomes Z-independent for even 
lower Z. 

Whereas M up '(Z) does depend on the specific stellar model code, the choice of physics 
used, and numerical implementation, as is well known, we expect that the general shape of 
the function arises from the underlying stellar physics and should not depend on the specific 
simulations used. 



4.3. Impact on supernova rates 



Using the ISalpeterl (119551 ) IMF we can compute the number of stars that make SNe, 



A/up' 



M" 2 ' 35 dM 
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where £o is a normalization constant. The lower limit, M up ', of the integral was explored 
here as a function of Z . A conservative upper limit for the maximum initial mass of a 
supernova progenitor, M max , may be as high as 100 M Q . Stars of higher mass are very 
rare. More reali stic is a typical transition mass from SNe to collapse without SN display at 



2011 



20 M - 30 M flHeger et al.ll2003l : lO'Connor & Ott 
16 M based on recent observational studies (jSmartt 

typical M max ps 20 M is not unreasonable, and we will use this as our reference case. Similar 



Ugliano et al.l 120121 ) to as low as 



20091 ) . Therefore a rough estimate of a 



to M up , there may not be a unique transition mass (IHeger et al.ll2003l ; lO ' Connor fc Ottl 



2011 



Ugliano et al.ll2012l ). Whereas the upper limit may depend on metallicity (e. g.. IZhang et al. 
20081) . this dependence is very uncertain and beyond the purpose of this paper. Woosley et al. 



(120021 ). Fig. 4, show that between 16 M and 27 M there is little difference in core mass, 
and hence anticipated SN properties, for stars of \Z\ = —4 and [Z] = and the same initial 
mass. Therefore we will consider here only a constant upper mass limit. 

Although the IMF itself could d epend on meta llicity as well, there is no overwhelming 
observational indication that it does (lKroupall2002l ). From theoretical considerations it was 
suspected tha t below a critical meta llicity Z CT ~ 10~ 6 only very massive (30 M — 300 M ) 
stars formed ( Ciardi fc Ferrara]|2005l ) that would not form supernovae and collapse to black 
holes instead. Only above Z a would a "nor mal" IMF set in. The recently reported very 
metal-poor (Z < 7xl0 -7 ) and low-mass star fjCaffau et al.ll201ll ) however calls into question 
the y alidity of Z CI . Recent work on nucleosynthesis from Pop III stars (IHeger fc Woosley 
20101 ) also does not provide strong support for a change in the IMF, at least not for massive 
stars. Hence we assume a metallicity-independent IMF. 

The Z-dependence of M up ' leads to a large change in Asn due to the high weight of 
the IMF at low M. Using our preferred value of M max = 20 M , our value M up ' = 8.35 M 
at [Z] = —2.3 (Table [Q this leads to a 25% increase of A SN ([Z] = —2.3) compared to 
A SN ([Z] = 0) with M up (Z Q ) = 9.2 M (Tabled]). Assuming the very conservative value 
of M max = 100 M for comparison, the change in Asn still is 20 %. The magnitude of 
the enhancement is not very sensitive to the upper mass limit as long as it is > 20 M . 
The uncertainty in M max dominates over the impact of the numerical noise in M up on the 
estimate of A SN . 

More important for cosmological applications is to know Asn as a function of redshift 
(ILien fc Fields! 120091 ). This would allow us to find the number of supernovae formed at 
each epoch of the history of the universe. To make this connection, we need to know 
the distribution of Z and star formation rate as a function of redshift. The s e para meters 
are all known to be connected to some extent; for example, iMannucci et al.l (120101 ) give a 
fundamental metallicity relation (FMR) that connects star formation rate, metallicity, and 
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galaxy mass, but no better than ±12 %. A detailed assessment of this, however, is beyond 
the scope of this work. The effect we present here is in co ntrast to the recently found lack 
of core collapse SNe at high redshift (jHoriuchi et al.ll201ll ). 



Similar to the SN rates, the variation of the lower mass limit with metallicity changes the 
ratio of neut ron star (NS) to black hole (BH) formation used in population synthesis studies 
(IFryerl Il999l ) . The effective nucleosynthesis yields from supernova used in galacto-chemical 
evolution models will be changed by allowing a larger contribution from low-met allicity stars. 



5. Conclusions 

We have computed for the first time the minimum mass for stars to make classical core 
collapse supernovae, M up ', as a function of metallicity, Z, at a high l e vel of detail. Our re sults 



generally agree with previous estimates by lEldridge Sz Toutl (120041 ); iHeger et al.l (120031 ) but 



significantly improve the quantitative results and detailed Z dependence. 

We find an increase of M up ' for [Z\ > —2, resulting in a decrease in the fraction of stars 
making CCSNe for a fixed IMF. Compared to a constant lower mass limit for supernovae 
this constitutes an increase in the supernova rate by 20 % to 25 % at low metallicity. Only 
for [Z] < —2 does the SN rate increase again. This has significant impact on interpreting 
observed SN rates as a probe of star formation history, on predictions of expected SN rates 
for upcoming surveys, and for NS/BH ratios for population synthesis studies. 

Future work should include stellar rotation, binary stars, the mass range required for 
ECSN from super-AGB stars, including how M up changes relative to M up ' as a function of 
metallicity, and a better understanding of the upper mass limit for supernovae. We need to 
connect the results to galactic parameters in order to better quantify the relative number of 
supernovae, neutron stars, and black holes as a function of redshift. Much of the underlying 
theory for these, however, is still very uncertain at the present. 
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Table 1. Summary of Model Fates 



z 




M u p' 


Z 




M up ' 


Z 




M u p' 






(M ) 






(M ) 






(M ) 







9.2 


1x10' 


-7 


8.8 


lxl0~ 3 




8.45 


1x10' 


-ii 


9.25 


1x10' 


-6 


8.7 


2.5x10- 
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8.65 


1x10' 


-10 


9.25 


4x10' 


-6 


8.56 


5xl0~ 3 




8.95 


3x10' 


-10 


9.35 


7x10' 


-6 


8.55 


7.5x10" 
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9.05 


5x10' 


-10 


9.5 


1x10' 


-5 


8.47 


1.5x10" 
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9.35 


1x10' 


-9 


9.5 


2x10' 


-5 


8.43 


1.9x10" 


2 


9.5 


3x10' 


-9 


9.3 


5x10' 


-5 


8.47 


3xl0~ 2 




9.6 


1x10' 


-8 


9.05 


1x10' 


-1 


8.35 


4xl0~ 2 




9.85 



